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摘 Æ: 以 2000 一 2019 年 SPOT/VEGETATION 的 NDVI 时 间 序 列 数 据 集 为 数据 源 , 采 用 年 均值 法 .Theil-Sen Medi- 
an 倾 斜 度 分 析 和 Mann-Kendall 检 验方 法 ,研究 宁夏 草地 植被 覆 被 的 时 空 分 布 及 变化 特征 ,利用 Hurst 指 数 方法 分 析 


草地 植被 覆 被 的 可 持续 性 特征 和 未 来 发 展 趋势 。 同 时 ,基于 地 理 探测 需 量 化 了 平均 降水 量 ,海拔 .地 区 生产 总 值 


(GDP) 等 13 个 因子 对 其 时 空 分 布 的 影响 。 结 果 表 明 :(1) 


2000—2019 年 间 ,宁夏 草地 植被 年 均 NDVI 呈 波动 增加 趋 


势 , 增 速 为 0.005'a ,区 域 波动 差异 性 较 大 ;空间 上 呈现 南 高 北 低 的 分 布 特征 , 极 高 植被 覆 被 和 高 植被 覆 被 区 域 集 中 


分 布 在 固原 市 六 盘山 地 区 以 及 沿 黄河 灌溉 地 带 ;(2) 20 a 间 植被 覆 被 状况 显著 改善 ,总 体 变化 趋势 向 好 ,但 仍 有 
59.341% 的 草地 未 来 可 能 存在 持续 退化 或 由 改善 向 退化 转变 的 潜在 风险 ;(3) 草地 植被 分 布 响应 最 敏感 的 环境 因子 
是 降水 ,交互 解释 力 整体 最 强 的 是 气候 与 土壤 ;影响 草地 植被 分 布 和 变化 特征 的 因子 之 间 交 互 作用 的 主要 表现 形 


式 为 相互 增强 或 非 线性 增强 关系 ,因子 之 间 不 存在 独立 关系 。 
关键 词 :NDVI; 趋势 分 析 ; 时 空 变化 ; 驱动 力 ; 草地 ; 宁夏 


草地 是 介 于 沙漠 与 森林 生态 系统 中 间 以 草地 
植被 占 主 体 的 土地 类 型 ,为 全 球 覆 盖 最 广 的 植被 类 
型 之 一 , 约 占 陆地 生态 系统 总 面积 的 1/4 一 。 草 地 
植被 在 保护 物种 多 度 .维持 畜牧 动态 平衡 .调节 水 
文 循环 和 气候 变化 等 方面 作用 重大 一 。 

归 一 化 植被 指数 (Normal Difference Vegetation 
Index, NDVD 与 地 表 植 被 叶 面 积 指 数 (Leaf Area In- 
dex, LAI) 净 初 级 生产 力 (Net Primary Production, 
NPP) 和 光合 有 效 辐 射 吸收 系数 (Fraction of Ab- 
sorbed Photosynthetically Active Radiation, FAPAR) 
等 密切 相关 ,可 用 来 准确 快速 评估 地 表 植 物 长 势 或 
HET, AR 地 表 及 人 类 活动 因素 对 草地 植被 影 
响 极其 显著 ,干旱 半 干 旱地 区 草地 植被 一 旦 遭 到 破 
坏 便 难以 恢复 “”。 因 此 ,研究 区 域 尺 度 的 草地 植被 
变化 及 其 驱动 因子 对 草地 生态 系统 的 保护 修复 具 
有 重要 意义 。 李 价 舌 等 ”对 影响 河套 灌区 林 草 植 
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被 覆盖 度 因素 进行 分 析 后 发 现 , 区 域 含水 量 的 大 小 
是 影响 河套 灌区 林 草 地 生长 状况 的 主要 驱动 力 ; 林 
妍 敏 等 ”发 现 ,宁夏 贺兰山 植被 覆盖 度 变 化 驱动 因 
子 主要 是 降水 .土壤 类 型 .植被 类 型 和 地 区 生产 总 
值 (CDP); 李 梦 华 等 ”研究 显示 ,气候 是 影响 宁夏 植 
被 覆盖 度 的 主导 因素 ; 付 含 培 等 ”对 黄河 流域 ND- 
VI 时 空 变化 及 驱动 力 进 行 分 析 得 出 20 a 来 黄河 流 
域 生态 环境 逐步 改善 ,气候 因素 是 影响 NDVI 空 间 
分 布 的 主导 因素 。 

目前 ,宁夏 植被 覆 被 变化 研究 中 的 植被 覆 被 包 
括 林地 .草地 湿地 耕地 .园地 和 其 他 土地 等 地 类 ， 
因 草 地 较 难 分 离 出 来 ,缺少 只 针对 草地 植被 覆 被 时 
空 动态 变化 及 驱动 力 的 研究 。 本 研究 采用 年 均值 
法 .Theil-Sen Median 倾斜 度 分 析 、.Mann-Kendall 检验 
方法 及 Hurst 指数 研究 宁夏 草地 植被 禾 被 时 空 分 布 
及 变化 特征 .可 持续 性 特征 和 未 来 发 展 趋势 ,同时 ， 
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应 用 地 理 探 测 需 (The geographical detectors ) 揭 示 影 
响 NDVI 变 化 的 驱动 力 。 地 理 探测 器 是 为 了 研究 空 
间 差 异 程度 及 其 显著 性 检验 和 归 因 分 析 , 广 泛 应 用 
于 医学 .考古 .土地 利用 、 区 域 经 济 等 各 个 领域 各 区 
域 尺 度 中 ,已 有 学 者 用 其 来 揭示 影响 植被 NDVI 变 
化 的 驱动 力 下 。 鉴 于 此 ,本 研究 利用 宁夏 2000 一 
2019 年 SPOTVEGETATION NDVI 时 间 序 列 数据 ， 
分 析 宁 夏 草地 NDVI 变化 特征 ,并 采用 地 理 探测 器 
进一步 分 析 宁 夏 草 地 植被 覆 被 变化 驱动 因子 ,为 宁 
夏 制定 科学 草地 保护 修复 模式 和 治理 措施 ,实现 草 
地 生态 系统 的 可 持续 经 营 提 供 科学 参考 和 决策 依 
据 , 对 于 维持 区 域 草 地 生态 系统 平衡 ,推动 黄河 流 
域 生态 保护 和 高 质量 发 展 具 有 重要 意义 。 


1 研究 区 概况 


宁夏 位 于 我 国 西北 地 区 (35° 14’ ~39。23' N, 
104°17’~107°39'E) (Al 1) , P5 i JE H AR a h 
腾 格 里 沙漠 . 乌 兰 布 和 沙漠 和 毛乌素 沙 地 相 围 , 南 
面 与 黄土 高 原 相连 。 地 势 南 高 北 低 , 地 貌 复杂 ,最 
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图 1 研究 区 示意 图 
Fig. 1 Study area schematic 
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高 点 为 敖包 疙 将 ( 俄 博 疙 郊 、 马 蹄 岭 ) ,海拔 3556 mo 
全 区 属 温带 大 陆 性 干旱 . 半 干 旱 气候 ,有 干旱 少雨 、 
蒸发 强烈 `. 风 大 沙 多 . 南 凉 北上 暖和 南 湿 北 干 等 特点 ， 
年 均 温 5~10 %C ,无 霜 期 150 d 左 右 ,年 降水 量 200~ 
700 mm, 是 我 国 典型 的 生态 脆弱 区 "*。 自 然 植被 类 
型 主要 有 草原 、 森 林 和 灌 从 等 。 宁 夏 草原 随 南北 气 
候 和 水 热 条 件 的 递 变 ,形成 了 丰富 多 样 的 草地 类 
型 ,由 南 向 北 依次 划分 为 :山地 草 多 Yk AE RE fa) 
原 、 温 性 典型 草原 . 温 性 荒漠 草原 `. 温 性 荒漠 Ti HE 
草原 化 荒漠 。2021 年 宁夏 林 草 资源 监测 结果 显示 ， 
温 性 荒漠 草原 占 草地 总 面积 的 62.113% , 温 性 典型 
草原 占 25.193% ,是 宁夏 草地 植被 面积 占 比 最 大 的 
两 个 草地 类 型 。 草 地 作为 宁夏 生态 系统 的 主体 ,发 
挥 着 防风 固沙 .保持 水 土 等 生态 功能 ,对 维持 区 域 
生态 平衡 ,促进 经 济 发 展 具有 重要 作用 。 


2 数据 与 方法 


2.1 数据 来 源 及 处 理 

人 研究 采用 的 数据 来 源 及 预 处 理 方 法 如 表 1 所 
示 。 在 建立 地 理 探测 器 模型 时 , 需 对 各 因子 数据 进 
行 预 处 理 , 其 中 地 貌 .土壤 数据 根据 一 级 类 别 ( 分 别 
为 地 貌 类 型 和 土 纲 ) 进 行 离散 化 处 理 ; 其 他 因子 均 
归 一 化 处 理 。 以 NDVI 为 Y 谈 量 ,13 个 因子 数据 为 
变量 进行 探测 。 
2.2 研究 方法 
2.2.1 时 间 序 列 与 空间 格局 分 析 逐年 计算 2000 一 
2019 年 宁夏 草地 NDVI 年 均值 ,分 析 20 a 间 NDVI 变 
化 规律 ; 逐 像 元 计算 2000 一 2019 年 草地 NDVI 平 均 
值 ,得 出 20a NDVI 平 均值 空间 分 布 特点 ,分 析 其 变 
化 特征 。 
22.2 _ NDVI 空间 波动 性 分 析 利用 变异 系数 反映 
草地 植被 覆 被 变化 的 波动 程度 ,计算 公式 "如 下 : 

1 1 < Ni 
C= ai = T2 (NDVI, - NDVI) (1) 
式 中 :C, 为 变异 系数 。 
2.2.3 Theil-Sen Median 倾斜 度 分 析 及 Mann-Kendall 
趋势 检验 ”Theil-Sen Median 倾斜 度 表 示 NDVI 的 变 
化 率 和 变化 趋势 “”。 计 算 公式 如 下 : 
NDVI,- NDVI, 
j-i 

Xt NDVI BY lal FEF (NDVL},¢ = 1,2,3, ,n, 

定义 斜率 5。 


S= Median ,1<i<j<n (2) 
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表 1 数据 来 源 与 处 理 


Tab. 1 Data source 


and data preprocessing 


数据 集 数据 来 源 空间 分 辨 率 数据 预 处 理 
NDVI 数 据 集 中 国 科学 院 资源 环境 科学 数据 1 km 一 
中 国 气 象 背景 数据 集 中 心 (https:Wwww.resdc.cn/) 500 m _ 
中 国 1:1000000 地 貌 类 型 空间 分 布 数据 1 km 一 
中 国土 壤 类 型 空间 分 布 数据 1 km 一 
(数字 高 程 模型 )DEM 90 m Al FA ArcGIS 10.8 获取 坡度 和 坡 向 数据 , 统 
一 范围 大 小 和 空间 分 辨 率 
GDP 密 度 宁夏 数据 网 站 (http://nxdata.com.cn/) HWE FE ArcGIS 10.8 中 ,插值 到 县 域 尺 度 
AORE 
草地 矢量 数据 2021 年 宁夏 林 草 资源 数据 库 一 按 细 化 地 类 提取 宁夏 回族 自治 区 范围 内 的 


式 中 :Median 为 中 位 数 ;NDVJL 和 和 NDVI 分 别 为 第 六 i 
时 间 序 列 的 NDVI 值 ; 当 斜 紊 S>0 时 ,NDVI 时 间 序 列 
{NDV 了 ) 呈 增长 趋势 ,反之 ,时 间 序 列 {NDVI)} 呈 下 降 
趋势 。 

Mann-Kendall 检验 (M-K 检 验 ) 在 评估 水 文 及 气 
候 变 化 趋势 中 广泛 应 用 , 它 考 虑 了 数据 的 有 效 分 
布 , 且 可 以 自动 处 理 异常 值 ,本 人 研究 用 来 判断 像 元 
尺度 上 NDVI 变 化 趋势 的 显著 性 。 

M-K 检 验方 法 计算 公式 参考 Ali 等 的 研究 。 
2.2.4 NDVI 未 来 演变 趋势 分 析 基于 Hurst 指 数 
法 ”对 宁夏 草地 NDVI 未 来 演变 趋势 进行 预测 ,其 

ANDVI,=NDVI,-NDVI,_,, i=1, 2,3,--,n (3) 


ANDVI,=1S.ANDVI, ,i=1, 2.3.2 (4) 
Lizi 


X(h)= X (ANDVI,- ANDVI), 1Shs: (5) 


i=l 
R(t)= max X(h)- min X(h), t= 1, 2,3,---,n (6) 


1 
t 


2 
S@= Savoy = ANDT}| t= 1,2,3,,n 


i=l 


(7) 
ROIS) = +t" (8) 
Th ce HAt AATE; H A Hurst H% WEA 
式 (8) 时 ,说 明 NDVI 具 有 长 期 记忆 的 时 间 序 列 , 即 
存在 Hurst 现象 ( 当 0.5 <H <1, RK AEE AE 
过 去 的 变化 )。 
2.2.5 草地 NDVI 趋 势 变化 驱动 力 分 析 选取 表 2 
因子 ,利用 地 理 探测 器 模型 分 析 影 响 宁 夏 草地 ND- 
VI 变化 潜在 驱动 力 。 采 用 ArcGIS 10.8 软件 将 影响 
因素 与 空间 位 置 进行 关联 ,分 析 NDVI 变 化 情况 。 


A 
草地 图 斑 ,融合 生成 草地 范围 矢量 数据 


表 2 宁夏 草地 NDVI 影 响 因素 
Tab.2 Influencing factors of grassland NDVI in Ningxia 


气候 因素 岂 表 因素 人 类 活动 因素 
平均 降水 量 (X,) EK (Xs) GDP(X,,) 
湿润 指数 (X;) 土壤 类 型 (X;) 人 口 密度 (Xs) 
“PIR BE (Xs) HUSK (Xo) 
年 均 温 (X,) 坡度 (Xo) 
20 CHA it (Xs) WE (Xa) 
210 CR (Xs) 


(1) 分 异 及 因子 探测 (g) 


SYNo? 
E (9) 
式 中 :为 YINDVIT) .X( 影 响 因 子 ) 的 分 类 ;h=1， 
2,°°,L; oy 为 h 层 方差 ; o 为 全 区 Y(NDVD) 的 方 
差 ; N, 为 h 层 单元 数 日 ;N 为 Y(NDVI) 单 元 数 日 。g 
的 值 域 为 [0,1], 根 据 g 值 大 小 可 以 验证 全 对 了 的 解 
释 性 ,解释 性 越 高 , 则 其 值 越 接 近 1。 

(2) 交互 作用 探测 

交互 作用 类 型 及 判断 准则 见 表 3。 


0 = 


表 3 探测 因子 交互 作用 类 型 


Tab.3 Interaction types of the detection factors 
交互 作用 类 型 判断 准则 
非 线性 减弱 qX NX) < Min[q(X)), q(X:)] 
单 因 子 非 线性 减弱 。 Min[g(X), gq(X3)] < qXNX) < Max[q(X,), qX)] 
双 因 子 增强 qX NX) > Max[q(X,), (OO 
非 线 性 增强 g(XiNX2) > g(Xi)+q(X2) 
独立 df IJ=4(CD)+dCD) 


(3) 生态 探测 () 


= NalNe = )NSSW,, 
J NaN r ~ SSW > 10) 
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L1 = k . oes Aa 
SSW, = ÈN, o? (11) 最 低 ,为 0.162;2018 年 NDVI 值 最 高 ,为 0.301。 宁 


h=1 
L2 

SSW) = > Nan (12) 
h=1 


式 中 :L1,L2 分 别 为 XY1,X2 分 层 的 数量 ; Na, Ny 分 
别 为 二 者 样本 数目 ; SSW,, , SSW, 分 别 为 Xl 和 X2 
形成 的 分 层 的 层 内 方差 之 和 。 

(4) 风险 区 探测 

用 于 探测 各 类 因子 对 植被 覆 被 影响 的 类 型 或 
适宜 范围 ,采用 :统计 量 检验 。 具 体 公 式 如 下 : 


et eae VY (13) 
RE 
Ny) Maer 


WP: Y, 表示 子 区 域内 的 NDVI 均 值 ; n, 为 子 区 域 

内 样本 数量 ;Var 表 示 方 差 。 统 计量 1 近似 地 服从 

Student” s 1 分 布 ,其 中 自由 度 的 计算 方法 为 : 
Var(Y,-,) i Var( _,) 


Nia Ni=2 


= 2 a 2 
1 Pen, 1 [va 
Myer ni npa l| ns 


(14) 
零 假设 成 : Y, 5Y, ,如 果 在 置信 水 平 a 下 拒 


H Ho, MAHAA F K ERE NDVI IHAT AE A H i 
的 差异 。 


3 结果 与 分 析 


3.1 草地 NDVI 时 间 序 列 变化 特征 

以 年 为 时 间 尺 度 ,分 析 宁 夏 草 地 NDVI 均 值 变 
化 情况 。 由 图 2 可 知 ,2000 一 2019 年 宁夏 NDVI 均 
值 呈 增加 趋势 , 增 速 为 0.005.a:。2000 年 NDVI 值 


df= 
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图 2 2000 一 2019 年 宁夏 草地 NDVI 年 际 变化 
Fig. 2 Interannual variation of grassland NDVI in Ningxia 
from 2000 to 2019 


夏 2000 一 2005 年 NDVI 均 值 为 0.199,2005 一 2010 年 
NDVI 44 {fi 4 0.214, 2010—2015 年 NDVI 均值 为 
0.253 , 2015—2019 年 NDVI 均 值 为 0.271, 草地 植被 
履 被 总 体 趋 势 向 好 发 展 。 
3.2 草地 NDVI 空 间 分 布 格局 

J&F 2000—2019 年 NDVI 时 间 序 列 数据 , 逐 像 
元 计算 20 a 平 均值 ,得 到 各 像 元 NDVI 平 均值 ,结合 
宁夏 实际 情况 ,本 人 研究 利用 自然 间断 点 分 级 法 (Jen- 
ks) ,将 草地 植被 状况 划分 为 5 种 类 型 : 极 低 植 被 覆 
#% (0 < NDVISO.190) 、 低 植被 覆 被 (0.190 < NDVIS 
0.251) .中 植被 覆 被 (0.251 < NDVI<0.338) .高 植被 
#8 BE (0.338 < NDVI<0.459) 、 极 高 植被 覆 被 (ND- 
VI>0.459)。 宁 夏 NDVI 空 间 分 布 总 体 表现 为 南 高 
北 低 ( 图 3)。 其 中 , 极 低 植被 覆 被 区 域 面积 ， 
7071.725 km’, 占 宁夏 草地 总 面积 的 34.397% ,主要 
分 布 于 宁夏 的 西北 部 地 区 ,地 处 腾 格 里 沙漠 边缘 ， 
包括 惠 农 区 ,大武 口 .中 宁县 .青铜 峡 、 沙 坡 头 ; 低 植 
被 覆 被 区 域 面积 为 7845.092 km?, H 38.157% ,多 位 
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图 3 2000 一 2019 年 草地 植被 覆 被 空间 分 布 
Fig. 3 Spatial distribution of grassland vegetation cover in 
Ningxia from 2000 to 2019 
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于 中 东部 干旱 少 十 地 区 ,包括 灵 武 市 .盐池 县 .利通 
区 \ 平 罗 县 、 兴 庆 区 、 永 宁县、 同心 县 、 红 寺 保 ;中 植 
被 覆 被 区 域 面积 为 3355.033 km’, 占 16.318% ,集中 
分 布 在 海原 县 、 贺 兰 县 ; 极 高 植被 覆 被 区 域 面积 》 
1688.375 km’, 占 8.212% , 高 植被 覆 被 面积 ; 
599.619 km’, 占 2.916% ,集中 分 布 在 固原 市 六 盘山 
地 区 以 及 沿 黄河 灌溉 地 带 。 
3.3 草地 NDVI 空 间 波动 性 分 析 

逐 像 元 计算 2000—2019 年 宁夏 草地 NDVI 变 异 
ARC) ,并 参照 自然 间断 法 进行 划分 (图 4)。 极 
低 波动 (0.041 < Cs<0.161) ,面积 1731.813 km’, 占 宁 
夏 草 地 面积 的 8.423% ; 低 波 动 (0.161 < C.<0.204) , 
面积 5501.385 km’, 占 26.758% ; 中 等 波动 (0.204 < 
C<0.241) , 面积 7538.200 km’, 占 36.665% ;高 波动 
(0.241 <C,<0.289) ,面积 4722.352 km’, 占 22.969%; 
极 高 波动 (C,>0.289) , 面积 1066.094 km’, 占 
5.185%。20 a 间 , 宁 夏 草 地 NDVI 变 化 整体 呈现 低 
波动 和 中 等 波动 变化 态势 ,区 域 波 动 差异 性 较 大 。 
3.4 草地 植被 时 空格 局 演变 

采用 Theil-Sen Median 倾斜 度 分 析 与 M-K 趋 热 
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图 4 宁夏 草地 NDVI 空 间 波动 性 
Fig.4 Spatial fluctuation of NDVI in Ningxia grassland 


显著 性 检验 相 结合 ,开展 草地 NDVI 变 化 趋势 分 
析 。 根 据 Theil-Sen Median 倾斜 度 分 析 结 果 , FES > 
0、S=0 和 S<0, 将 草地 植被 覆 被 变化 分 别 划 分 为 改 
He .基本 不 变 和 退化 3 种 情况 。 根 据 M-K 检 验 结 
(0.05 置信 水 平 ) 将 NDVI 变化 显著 性 分 为 5 级 。 将 
Theil-Sen Median 倾斜 度 分 级 结果 与 M-K 检 验 分 级 
结果 进行 空间 闭 加 ,将 草地 NDVI 变 化 趋势 划分 为 : 
显著 改善 ,轻微 改善 .基本 不 变 、 轻 微 退化 和 显著 退 
化 5 类 中 。 由 变化 趋势 结果 可 见 , 宁 夏 草地 20 a 间 
植被 覆 被 状况 明显 改善 ,退化 面积 较 小 ,总 体 变 化 
趋势 向 好 。 其 中 ,显著 改善 面积 13734.576 km’, 占 
66.803% ;轻微 改善 面积 5116.118 km’, rF 24.884% ; 
基本 不 变 面 积 84.041 km?, 占 0.409% ;轻微 退化 面积 
1575.062 km , 占 7.661%; 显著 退化 面积 50.047 km’, 
占 0.243%( 图 5)。 
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5 2000 一 2019 年 宁夏 草地 植被 覆 被 变化 趋势 
Fig. 5 Change trend of grass coverage in Ningxia 
from 2000 to 2019 


为 进一步 分 析 宁 夏 草 地 演化 过 程 ,利用 ArcGIS 
10.8 进行 NDVI 空 间 分 析 ,得 到 2000 一 2019 年 宁夏 
草地 植被 覆 被 空间 转移 矩阵 ( 表 4) ,结果 表明 ， 
2000—2019 年 期 间 ,草地 极 低 植被 覆 被 大 幅 减 少 ， 
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表 4 宁夏 草地 植被 覆 被 空间 转移 矩阵 


Tab. 4 Spatial transfer matrix of grassland vegetation cover in Ningxia I% 
2000 4 
植被 类 型 i a TEE EEEE es 5 
极 低 植被 覆 被 低 植被 覆 被 中 植被 覆 被 高 植被 覆 被 极 高 植被 覆 被 合计 

2019 年 极 低 植被 覆 被 15.746 0.040 0.007 0.010 15.803 
低 植 被 覆 被 29.591 0.297 0.008 29.896 

中 植被 覆 被 28.450 1.291 0.071 0.022 29.834 

高 植被 覆 被 9.923 4.493 0.692 0.086 0.021 15.215 

极 高 植被 覆 被 0.304 3.093 3.477 1.735 0.643 9.252 

合计 84.014 9.214 4.255 1.831 0.686 100.000 

104°0'E 105°0'E 106°0'E 107°0'E 


低 一 级 植被 覆 被 向 高 一 级 植被 覆 被 转化 明显 。 总 
体 来 看 ,宁夏 草地 植被 增长 明显 ,生态 保护 修复 成 
效 显 著 。 

3.5 草地 NDVI 未 来 演化 趋势 

3.5.1 NDVI 空 间 可 持续 性 经 R/S 分 析 , 宁 夏 草地 
Hurst 指数 均值 为 0.488, 说 明 全 区 草地 植被 覆 被 变 
化 存在 时 间 序 列 上 反 持 续 性 。 分 类 结果 显示 ,0.5 < 
<1 的 区 域 面积 合计 8884.793 km’, 占 全 区 草地 总 
面积 的 43.214% ,该 区 域 植 被 履 被 呈 同 向 持续 性 特 
征 , 即 植被 覆 被 变化 可 能 维持 2000 一 2019 年 期 间 的 
变化 趋势 ;0< 且 <0.5 区 域 面积 合计 11675.051 km’, 
占 56.786%, 植 被覆 被 旦 反问 持续 性 特征 , 即 植被 覆 
被 未 来 变化 可 能 与 2000 一 2019 年 期 间 的 趋势 相反 。 
3.5.2 NDVI 未 来 发 展 趋 势 将 NDVI 空 间 变 化 趋势 
和 空间 可 持续 性 特征 进行 空间 麦 加 分 析 ( 图 6) , 根 
据 分 析 结 果 将 宁夏 草地 植被 覆 被 未 来 发 展 趋势 分 
为 良性 方向 恶性 方向 稳定 不 变 和 不 确定 4 个 方向 
( 表 5)。 未 来 草地 植被 覆 被 可 能 呈 良 性 发 展 的 区 域 
面积 占 草地 总 面积 的 40.255% , 主要 分 布 在 原 州 区 、 
PAE 泾 源 县 、 隆 德 县 西 吉 县 .盐池 县 东南 部 E] 
心 县 北部 海原 县 北部 .中 宁县 南部 ` 沙 坡 头 区 东南 
部 和 大 武 口 区 。 未 来 草地 植被 覆 被 可 能 呈 恶 性 发 
展 的 区 域 面积 占 草 地 总 面积 的 59.341% , 主要 分 布 
在 兴 庆 区 、 惠 农 区 、 永 宁县、 同心 县 南部 海原 县 南 
部 ` 沙 坡 头 区 北部 .中 宁县 北部 利通 区 .红寺堡 ER 
池 县 西北 部 、 灵 武市 ;未 来 草地 植被 履 被 发 展 方向 
稳定 不 变 和 不 确定 的 区 域 合计 面积 占 草 地 总 面积 
的 0.404%。 

3.6 草地 NDVI 空 间 分 布 的 影响 因子 

3.6.1 不 同 影 响 因子 对 NDVI 空 间 分 异 的 驱动 力 分 
析 采用 地 理 探测 器 对 NDVI 的 空间 分 异 及 变化 趋 
势 进行 分 析 ,因子 探测 结果 见 表 6, 不 同 因子 对 宁夏 
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注 :aa 为 持续 性 & 显 著 退 化 ; ab 为 持续 性 & 轻 微 退化 ; ac 为 持续 性 区 
稳定 不 变 ; ad 为 持续 性 区 轻微 改善 ; ae 为 持续 性 & 显著 改善 ; ba 为 
反 持 续 性 区 显著 退化 ; bb 为 反 持续 性 区 轻微 退化 ; be 为 反 持续 性 & 
稳定 不 变 ; bd 为 反 持续 性 皮 轻 微 改善 ; be 为 反 持 续 性 区 显著 改善 。 
图 6 宁夏 草地 植被 覆 被 未 来 发 展 趋势 


Fig.6 Future development trend of grass coverage in Ningxia 


草地 NDVI 空 间 分 布 的 解释 能 力也 不 同 。 其 中 ,年 
平均 降水 量 、 湿 润 度 .干燥 度 .年 均 温 .积温 (>=0%C ) 
的 解释 力 均 超过 60% ,分 别 为 80.0% .79.4% .77.0% 、 
62.0% .61.5% ,是 主要 影响 因子 ;海拔 的 解释 力 为 
59.1% ; 土壤 类 型 .地貌 .坡度 . 坡 向 、 积 温 (=10 C), 
AO GDP 对 宁夏 草地 NDVI 影 响 不 显著 , 均 无 解 
EH. 
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表 5 宁夏 草地 植被 覆 被 未 来 发 展 方向 面积 占 比 
Tab. 5 Distribution of future development direction of 
grass coverage in Ningxia 
发 展 方向 结 面积 /km 占 比 /% 
良性 方向 持续 性 & 轻 微 改 善 1554.343 7.560 
持续 性 & 显 著 改 善 5694.847 27.699 
反 持 续 性 & 轻 微 退 化 1009.142 4.909 
反 持 续 性 & 显 著 退 化 17.953 0.087 
恶性 方向 持续 性 & 轻 微 退 化 565.989 2.753 X, 
持续 性 & 显 著 退 化 32.126 0.161 Xo 
反 持 续 性 & 轻 微 改善 3562.232 17.326 © x, 
反 持 续 性 & 显 著 改 善 8039.117 39.101 Xu 00000 ° ee @ o v, 
不 确定 反 持 续 性 & 不 变 51.024 0.248 Xia eee © S OO «|x, 
稳定 持续 性 & 不 变 33.071 0.156 Xa (1 11 © e@ Se @ 
表 6 因子 探测 结果 0 0.1 0.2 0.29 0.39 0.49 0.58 0.68 0.78 0.88 1 
Tab. 6 Results of factor detector 
7 lil PE S E ee 用 由 强 到 弱 。XX 为 年 平均 降 
EERO) ca use 水 量 oe aa Xs J20 积温 ,Xo 为 
THEO) orod nas 10 CHA XAT TEAK , Xe J ERY, Xo THA, Xo ABE Xi 为 坡 
FEE) ITO dais H], Xa J GDP, X HA ORE. FEl 
EHX) 0.620 0.020 图 7 交互 作用 探测 影响 力 
=0 CHUL) 0.615 0.014 Fig. 7 Interactive detection explanatory power 
HHR) 0.591 0.031 of various factors 
土壤 类 型 (X;) 0.499 0.564 
SHA (Xo) 0.392 0.870 X% X 
AHERE Xs) 0.258 0.387 % NX 
YUE (Xio) 0.179 0.576 X N N 7 
H (Xu) 0.170 0.868 x, aa X, 
=10 CHX) 0.137 0.880 x Y Y N N 4s 
GDP(X::) 0.071 0.908 x: ees Xs 
xX, aaa ee ea 7 
3.6.2 交互 作用 探测 分 析 绝 大 多 数 因子 交互 呈现 。 Y Y Y NNNN % 
相互 增强 ,表现 为 两 因子 的 4 值 和 大 于 单 因 子 / 值 。 宇 YY YSN ANN 
pom OO OY YN OY ON i ES 
We A PPR AGAR AT A, FLA XX 130.975 )>XXs EY YY eo a 
(0.937)>XsUXs(0.935)>Xs1%s(0.930)>X1%(0.928) T 和 
交互 作用 影响 力 排 名 前 五 ,土壤 与 以 人 口 密度 为 代 
注 :置信 水 平 为 95%,Y 表 示 两 因子 对 NDVI 空 间 分 布 有 显著 差异 ， 


表 的 人 类 活动 因子 交互 作用 最 强 , 但 土壤 与 以 降 
水 .气温 为 代表 的 主导 因子 交互 作用 整体 表现 更 为 
明显 (图 7)。 虽 然 , 人 类 活动 的 单 因 子 解释 力 较 低 ， 
但 与 环境 因子 的 交互 作用 解释 力 较 强 ,表明 人 类 活 
动因 素 是 影响 植被 覆 被 空间 分 异 格局 的 间接 因素 。 
3.6.3 生态 探测 分 析 由 图 8 可 知 ,大 部 分 交互 因 
子 对 NDVI 空 间 分 布 有 显著 性 差异 ,年 平均 降水 量 
及 湿润 度 与 其 他 因子 间 更 为 突出 。 由 此 可 见 , NDVI 
受气 候 因 素 .地 表 因 素 及 人 类 活动 因素 共同 影响 ， 


N 则 相反 。 
图 8 生态 探测 结果 


Fig. 8 Ecological detection results 


气候 因素 为 主导 。 

3.6.4 风险 探测 分 析 ”基于 风险 探测 ,可 以 得 出 区 
域 草 地 NDVI 值 最 高 时 的 各 类 影响 因子 的 范围 ( 表 
7)。 气 候 因素 :年 平均 降水 量 545~623 mm ,年 均 温 
2.8~3.7 % ,表明 降雨 量 大 气温 低 的 地 区 草地 植被 
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表 7 不 同 影响 因子 的 适宜 类 型 或 范围 


Tab.7 Suitable ranges or types of detection indicators 


因子 草地 NDVI 适 宜 类 型 或 范围 NDVI 均 值 
年 平均 降水 量 545~623 mm 0.596 
湿润 度 80~200 0.596 
干燥 度 0~62 0.542 
年 均 温 2.8~3.7 C 0.472 
0 CA 2165~2416 C 0.472 
=10 % 积 温 1905~2087 C 0.402 
海拔 2303~2502 m 0.436 
土壤 类 型 半 水 成 土 0.403 
地 貌 中 起 伏 山 脉 0.596 
坡度 24.73°~28.73° 0.365 
坡 向 349.69°~0.84° 0.385 
GDP 5x10'~6x10°96 km 0.475 
人 口 密度 62~72 人 :km? 0.587 


履 被 状况 更 好 。 地 表 因 素 :高程 2303~2502 m ,坡度 
24°~28° ,土壤 类 型 为 半 水 成 土 ,地 貌 为 中 起 伏 山脉 
时 草地 植被 覆 被 状况 最 好 ,说 明 在 高 山地 区 草地 植 
被 生长 状况 更 好 。 人 类 活动 因素 :GDP(Sx10'~6x 
10' 元 km“) 、 人 口 密度 (62~72 人 .km2) 时 草地 植被 
履 被 状况 更 好 ,说 明 受 人 为 干扰 小 的 地 区 草地 植被 
履 被 状况 最 好 。 


4 讨论 


4.1 宁夏 草地 NDVI 时 空 变化 特征 

宁夏 位 于 中 国 西北 部 ,地 处 黄河 中 上 游 地 区 及 
腾 格 里 沙漠 .黄土 高 原 和 蒙古 高 原 的 交界 地 带 ,地 
形 复杂 ,气象 .地 表 和 人 类 活动 等 因素 共同 影响 区 
域内 草地 植被 覆 被 的 变化 。 从 时 间 上 看 ,2000 一 
2019 年 宁夏 草地 年 平均 NDVI 值 呈 增 加 趋势 , 增 速 
为 0.005.a, 宁 夏 2000 一 2019 年 NDVI 值 由 0.162 上 
升 至 0.290。21 世纪 初 ,国家 对 黄土 高 原 实 施 了 退 
耕 还 林 ( 草 ) 工 程 等 一 系列 生态 环境 治理 措施 ,有 效 
促进 了 植被 生长 和 恢复 。 宁 夏 在 此 期 间 全 域 坚 持 
天 然 草原 禁 牧 封 育 措施 ,植被 生长 良好 ,草原 植被 
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的 气候 则 不 利于 植被 的 生长 1。 

20a 间 ,宁夏 草地 NDVI 空 间 分 布 变化 整体 呈现 
低 波 动 和 中 等 波动 变化 态势 ,区 域 波 动 差异 性 较 
大 。 低 波动 和 极 低 波动 区 域 主要 分 布 于 南部 山区 
降水 相对 充沛 区 域 和 黄河 沿岸 灌溉 便利 地 区 ;中 高 
波动 区 域 主要 分 布 于 中 部 干旱 带 .荒漠 地 区 以 及 北 
部 贺兰山 地 区 等 ,这 与 宁夏 降水 量 年 度 波动 情况 密 
切 相 关 。 

目前 ,宁夏 草地 显著 退化 面积 仅 占 草 地 面积 的 
0.243% , 但 根据 未 来 发 展 趋势 分 析 结 果 仍 有 
59.341% 的 草地 未 来 可 能 存在 持续 退化 或 由 改善 向 
退化 转变 的 潜在 风险 。 林 妍 敏 等 "研究 发 现 ,宁夏 
贺兰山 地 区 仍 有 30.7% 的 植被 存在 退化 的 潜在 风 
今 。 草 地 植被 覆 被 对 气候 变化 的 响应 较 敏感 ,气候 
条 件 是 影响 该 区 域 草地 植被 生长 的 主要 因素 3, 也 
有 学 者 指出 干旱 半 干 旱 区域 , 植 被 生长 受气 候 因 素 
限制 , 暧 湿 化 有 利于 植被 生长 ”。 草 原 植被 恢复 主 
要 受气 象 条 件 和 生态 保护 政策 双重 因素 的 影响 ,未 
来 应 加 强 草 原 保 护 修复 工程 建设 的 长 期 性 和 连续 
性 ,建议 继续 推行 禁 牧 及 退耕 还 林 还 草 政策 ,建立 
健全 草地 保护 利用 相关 法 规 体 系 ,深化 牧区 改革 ， 
优化 畜 种 冀 群 结构 ,控制 草地 载 畜 量 。 草 原 管 理 部 
门 要 加 大 草地 管理 力度 ,进一步 加 强 对 农 牧 民 的 宣 
传教 育 , 杜 绝 乱 占 、 乱 暴 、 乱 牧 等 现象 。 同 时 ,要 完 
善 降水 .气温 等 气候 条 件 的 动态 变化 监测 工作 ,及 
时 掌握 气候 变化 对 草原 植被 生长 的 影响 ,为 草地 保 
护 和 培育 决策 提供 科学 支撑 。 
4.2 驱动 因子 对 草地 NDVI 空 间 分 布 的 响应 

草地 生态 系统 对 气候 因素 变化 响应 尤为 显著 ， 
尤其 是 在 半 干 旱地 区 ,气候 是 对 草地 植被 生长 具有 
决定 性 影响 的 自然 驱动 力 ? 六 。 本 研究 分 析 了 草地 
NDVI 与 气象 .地 表 、 人 类 活动 因素 之 间 的 关系 , 结 
果 表 明 ,年 平均 降水 量 .湿润 度 .干燥 度 .年 均 温 、 积 
温 (=0%) 的 解释 力 均 超过 60% ,其 中 年 平均 降水 
量 与 湿润 度 的 解释 力 均 接近 80%。 由 此 可 见 , 降 水 


群落 趋 于 稳定 下 ,草原 生态 得 到 明显 恢复 。 从 空 
间 上 看 ,宁夏 草地 NDVI 分 布 总 体 表 现 为 : 南 高 北 低 ， 
与 纬度 负 相 关 , 空 间 分 异 特征 明显 。 王 文 晓 等 *“ 研 
究 中 也 发 现 ,宁夏 草原 物种 丰富 度 由 南 到 北 呈 递减 
趋势 。 这 与 区 域 的 气温 与 降水 等 有 着 密切 关系 , 宁 
夏 具 有 南 寒 北 暖 . 南 湿 北 于 特点 ,气温 大 体 呈 北 高 
南 低 分 布 ,降水 量 由 南 向 北 递减 ,而 高 气温 低 湿 度 


是 宁夏 草地 植被 覆 被 空间 分 异 的 主要 驱动 力 , 这 与 
李 梦 华 等 ”研究 结果 一 致 。 本 研究 得 出 土壤 与 气 
象 因素 中 的 降水 .气温 交互 作用 表现 更 为 明显 , 原 
因 是 宁夏 土壤 有 机 质 自 南 向 北 呈 降低 趋势 , 南 壤 北 
沙 分 布 格局 明显 ”。 风 险 探测 分 析 结 果 显 示 ,土壤 
类 型 为 半 水 成 土 区 域 草地 植被 覆 被 最 好 。 半 水 成 
土 是 指 在 气候 比较 湿润 的 森林 草原 条 件 下 发 育 而 
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成 的 土壤 ,包括 草 旬 土 .黑土 潮 土 AKE WME 
等 中 , 相 比 干旱 土 更 适宜 植被 生长 。 刘 完 锋 等 " 指 
出 ,在 干旱 半 干 旱地 区 ,土壤 类 型 作为 植被 生长 和 
空间 分 布 的 重要 生态 因子 ,对 植被 生长 和 雨水 再 利 
用 效率 的 影响 显著 。 

本 研究 采用 地 理 探测 器 对 宁夏 草地 NDVI 空 间 
分 异 特征 进行 分 析 ,结论 与 同 区 域 研究 结果 基本 一 
致 *”™。 利 用 地 理 探测 器 探究 驱动 因子 优势 明显 ， 
可 以 从 空间 上 定量 研究 各 数值 型 与 定性 数据 驱动 
因子 对 研究 区 草地 NDVI 演 变 的 影响 ,但 人 类 活动 
对 于 草地 植被 覆 被 的 影响 十 分 复杂 ,特别 是 管理 模 
式 ` 家 畜 结 构 与 密度 .放牧 强度 等 均 会 对 草地 植被 
生长 造成 一 定 影响 。 同 时 ,草原 火灾 、 鼠 虫 灾害 等 
难以 量化 的 因子 也 会 对 草地 植被 覆 被 造成 影响 ,以 
上 因素 本 研究 尚未 考虑 。 建 议 今后 研究 中 可 以 选 
取 其 他 影响 草地 植被 覆 被 的 相关 因子 ,深入 探究 植 
被 覆 被 状况 变化 的 驱动 机 制 。 


5 结论 


本 研究 利用 2000 一 2019 年 SPOT/VEGETA- 
TION NDVI 数 据 、 气 象 数据 等 ,采用 年 均值 法 、Theil- 
Sen Median 倾斜 度 分 析 、Mann-Kendall 检 验方 法 和 
地 理 探测 器 模型 ,分 析 了 宁夏 草地 植被 覆 被 时 空 变 
化 特征 及 其 影响 因子 的 驱动 作用 。 结 论 如 下 : 

(1) 20 a 间 ,宁夏 草地 植被 覆 被 呈现 波动 上 升 
趋势 ,植被 覆 被 状况 整体 呈现 明显 改善 的 趋势 ; 宁 
夏 草地 植被 覆 被 在 空间 分 布 上 总 体 呈 现 南 高 北 低 
的 特征 ,高 覆 被 区 集中 在 固原 市 六 盘山 地 区 以 及 沿 
黄河 灌溉 地 带 ; 空 间 分 布 变化 整体 以 中 等 波动 和 低 
波动 变化 态势 为 主 , 区 域 波 动 差 异性 较 大 。 

(2) 宁夏 草地 植被 覆 被 在 2000—2019 年 间 虽 
明显 改善 ,但 预测 未 来 发 展 趋势 持续 性 较 差 , 仍 有 
59.341% 的 草地 未 来 可 能 存在 持续 退化 或 由 改善 向 
退化 转变 的 潜在 风险 。 

(3) 宁夏 草地 植被 履 被 的 空间 分 布 受 气象 .地 
表 和 人 类 活动 等 因素 共同 作用 ,气象 因素 影响 作用 
大 于 地 表 与 人 类 活动 因素 ,年 平均 降水 量 对 植被 覆 
被 空间 分 布 的 影响 最 大 。 

(4) 影响 草地 植被 分 布 和 变化 特征 的 因子 之 间 
交互 作用 的 主要 表现 形式 为 相互 增强 或 非 线性 增 
强 关 系 , 因 子 之 间 不 存在 独立 关系 ,交互 解释 力 最 
强 的 是 气候 与 土壤 。 
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Analysis of the spatiotemporal variation characteristics and driving factors of 
grassland vegetation cover in Ningxia based on geographical detectors 
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National Forestry and Grassland Administration on Ecological Hydrology and Disaster Prevention in Arid Regions, 
Xi’an 710048, Shaanxi, China) 


Abstract: This study aims to examine the spatiotemporal variation characteristics of grassland vegetation cover 
at the regional scale and analyze its driving factors. The findings will provide a scientific reference and decision- 
making basis for the scientific formulation of protection and restoration models, treatment measures, and the 
sustainable management of the grassland ecosystem in Ningxia, which are crucial for maintaining the balance of 
the regional grassland ecosystem and promoting ecological protection and high- quality development in the 
Yellow River Basin. In this study, the NDVI time series dataset of SPOT/VEGETATION (2000-2019) was used 
as the data source. The annual mean method, Theil- Sen Median trend analysis, and Mann- Kendall test were 
employed to study the spatiotemporal distribution and variation characteristics of grassland vegetation cover in 
Ningxia. Furthermore, the Hurst index method was used to analyze the sustainability characteristics and future 
development trends of grassland vegetation cover. Simultaneously, the influence of 13 factors, such as average 
precipitation, altitude, and gross domestic product, on the spatiotemporal distribution was quantified based on the 
geographical detectors approach. The results show that from 2000 to 2019, the average annual NDVI of 
vegetation in Ningxia grassland showed a fluctuating growth trend, with a growth rate of 0.005 per year. The 
regional fluctuation was quite different, with extremely high and high vegetation cover areas concentrated in the 
Liupan Mountains and the irrigation area along the Yellow River. Overall, the NDVI change showed a low to 
medium fluctuation trend, and the regional fluctuation was quite different. The vegetation cover condition 
improved significantly over the 20-year period, with a small degradation area and a favorable overall change 
trend. However, 59.341% of the grasslands are projected to face potential risks of continuous degradation or 
transformation from improvement to degradation in the future. The most sensitive environmental factor 
influencing grassland vegetation distribution response was precipitation, and climate and soil had the strongest 
interaction explanatory power overall. The relationship between the factors affecting the distribution and 
variation characteristics of grassland vegetation primarily manifested as mutual reinforcement or nonlinear 
enhancement, with no independent relationship between the factors. This study provides a scientific reference and 
decision-making basis for the sustainable management of the grassland ecosystem in Ningxia. 
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